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Abstract 

Macaulay dual spaces provide a local description of an affine scheme 
and give rise to computational machinery that is compatible with the 
methods of numerical algebraic geometry. We introduce eliminating dual 
spaces, use them for computing dual spaces of quotient ideals, and develop 
an algorithm for detection of embedded points on an algebraic curve. 


1 Introduction 

We explore algorithms that use dual spaces to compute local information about 
polynomial ideals over the complex numbers, as an alternative to working di¬ 
rectly with polynomials. This strategy often has advantages in the regime of 
approximate numerical computations and for hybrid symbolic-numerical algo¬ 
rithms. Many computations are reduced to linear algebra, allowing numerical 
linear algebra techniques to be applied. 

Given generators of a polynomial ideal / and a point p in its vanishing set, 
the dual space of I at p is the vector space dual of the extension of I in the 
local ring at p, and it uniquely encodes the local properties of / there. Certain 
combinatorial information about the dual space, such as dimension, can be 
accurately computed even when p is only known approximately but with high 
enough precision. Existing methods of numerical algebraic geometry (see e.g., 
PD] and m) offer efficient algorithms to approximately compute points on the 
vanishing set of an ideal. 

The contributions of this paper are: 

• Establishing a correspondence between a local monomial order on a local 
ring (primal order) and the corresponding dual order on the monomials of 
the dual space. 

• Characterizing the dual space of colon (quotient) ideals. 
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• Introducing the notion of eliminating dual spaces as a computational tool. 

• Applying these ideas to construct an algorithm for detecting embedded 
primary components of a curve. 

Additionally, we consolidate the necessary background on dual spaces, state 
their properties in a modern language, and provide proofs where references are 
unavailable. 

The idea of studying systems of polynomials through dual spaces dates back 
to Macaulay mi- Most of the recent work using Macaulay’s machinery con¬ 
cerns zero-dimensional ideals or, geometrically speaking, isolated points. This 
includes algorithms for computing a basis of the dual space ns is and the 
local Hilbert function at an isolated point [6] , as well as various deflation proce¬ 
dures mmm- Several studies depart from the zero-dimensional setting: the 
local dimension test [3, computations using dual spaces for homogeneous ide¬ 
als and an numerical algorithm for local Hilbert polynomial in the general 
case US¬ 
As an application of eliminating dual spaces we present an algorithm for 
detecting embedded components of a curve. This fills in a key missing piece 
in the program to numerically compute a primary decomposition of a polyno¬ 
mial ideal as laid out in na in the case of one-dimensional ideals. We remark 
that eliminating dual spaces may, in theory, be used to provide an embedded 
component test in dimension more than one. However, we did not pursue this 
direction since for practical computation this potential technique appears to be 
inferior to an alternative treatment we develop in |13j . Our algorithms for dual 
spaces and detecting embedded components are implemented in the Macaulay2 
computer algebra system [4] and the code is posted at [121 . 

The rest of the paper is structured as follows. Section [2] gives an expository 
background on dual spaces of ideals in a local ring, laying out the facts that will 
be used later in the paper and giving proofs where they were not readily available 
in the literature. Additionally we characterize the dual space of quotient ideals. 
Section [3 describes the correspondence between the initial terms of an ideal 
with respect to a (primal) local monomial order and the initial terms of the 
dual space in the dual order. Such a relation has been observed before in 
special settings such as zero-dimensional ideals with a graded order, but we give 
a fully general result. This section also describes the relation between truncated 
dual spaces and the Hilbert function. Section [I] introduces the eliminating dual 
space, which generalize the truncated dual space, and can be used to compute 
truncated dual spaces of a quotient ideal. These ideas are applied in [5] to the 
problem of detecting the embedded components on a curve, and examples are 
provided. 


2 Macaulay dual spaces 

For a £ (Z> 0 ) N and y £ C^, define 
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• Oi\ = QiSo^! • • • 

• d a = i2!^, and 

a! 7 

• 9“[y] : i? —> C defined by d a [y\(g) = ( d a g)(y). 

The differential functional d a [y\ sometimes would be written d xa [y\ (e.g. d 1 — 
d y + d x yz ) and when the point y is implied d a [y\ would be written as d a . For 
y G C N , let D y = span c { d a [y\ | a G (Z> 0 ) N ) be the vector space of differential 
functionals at y. This linear space is graded by order , for a finite sum q = 
J2c a d a , 

ord q = max led. 
c a ^ o 

The homogeneous part of order i of q G D y is referred to as qi. This grading is 
the associated graded linear space of the filtration D*: 

Dy C Dy C Dy C , where D '* = {q G D y | or d q < i}}. 

Definition 2.1. The Macaulay dual space, or simply dual space, of differential 
functionals that vanish at y for an ideal I C C[x] = C[xi,..., xjv] is 


D y [I] = {q G D y | q{g) = 0 for all g G I}. 


(1) 


The dual space D y [I] is a linear subspace of D y , a basis of D y [I] is called a dual 
basis for I. 

The following theorem of Macaulay describes the dimension of the dual space 
at an isolated solution y. The following statement appears in the classical text 
of Macaulay El- 

Theorem 2.2. A solution y G V(J) is isolated with multiplicity m if and only 
if dime D y [I} = to. 

Definition 2.3. A subspace S C D y is homogeneous if it is spanned by ho¬ 
mogeneous functionals q G D° rAq \ Z?° rd9_1 . If, in addition, S is spanned by 
d a \y\, a G A for some subset A C (7j>q) n , then S is called monomial. 


2.1 Local ring vs. Dual space 

For the purpose of this section, without a loss of generality, we may assume 
y = 0 G C N . Consider the local ring R 0 = R m where m = {x\, ..., Xi v). Let the 
space of dual functionals be defined as above replacing R (polynomial) with Rq 
(rational functions with denominators not vanishing at 0). 

Remark 2.4. Ideals in R with no primary components away from the origin 
are in one-to-one correspondence with ideals in the local ring Rq: 

• an ideal I C R defines the extension IRq C Rq; 
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• an ideal I C Rq corresponds to the ideal I fl R C R with no primary 
components away from the origin. 

Proposition 2.5. For ideal I C R, the dual space Dq[I\ is identical to the dual 
space of its extension in Rq, Dq[IRq], 

Proof. Any rational function g £ Rq can be expressed as a power series g = 
^ a c a x a . If q £ Do[I\, then for any / £ /, q(x a f) = 0 for all monomials x a . 
Then 

9 ( 5 /) = c a.q{x a f) = 0 

<y.£(K> 0 ) N 

so q £ D 0 [IRo\. For q not in D 0 [I\ there is some f € I with q(f) ^ 0, and / is 
also in IRq. □ 

As a result we will speak interchangeably about the dual space of an ideal I 
at the point 0 and the dual space of its extension in the localization of R at 0, 
IR q . 

The following lemma provides another characterization of the extension of 
an ideal I in the local ring, which will help describe the close connection between 
IRq and the Macaulay dual space. 

Lemma 2.6. For any ideal I C R, 

OO 

IR 0 nR = P|(/ + m fc ). 

k =1 

Proof. Let R denote the completion of R with respect to the maximal ideal 
m (the formal power series ring R = C[[xi,... , xjv]])- The kernel of the map 
of .R-modules R/I R/I is P| fc m k (R/I), and by the exactness of completion 
R/I = R/IR (see [I] Chapter 10). Composing the quotient map R —> R/I with 
the above, we see that IR fl R, which is the kernel of natural map R —> R/IR , 

is n fc 0r + m fe ). 

For any / £ IR fl R, there is h € /, g £ R such that / = hg , so g = h/f is a 
rational function in R 0 . Therefore / £ IRq fl R, and so IRq fl R = IR flf?. □ 

Proposition 2.7. For ideal I C R, f £ IRq fl R if and only if q(f) = 0 for all 
q £ Dq[I}. 

Proof. It follows from the definition that / £ I implies q(f) = 0 for all q £ Dq[I\. 

Let R k be the space of polynomials with degree < k, let f k denote the 
truncation of / to degree k and let I k C R k be the set {f k : f £ I}. Since R k 
is a finite dimensional vector space, 

(/ fe ) x =D k [I\ = D 0 [I + m k+1 }. 

Suppose for some polynomial / that q(f) = 0 for all q £ D 0 [I] = 

Because ( I k )- L± = I k , we have f k £ I k which implies f £ I + m fc+1 . By 
Lemma [2~dl / £ IR 0 fl R. □ 
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Corollary 2.8. For ideals J\, J 2 C Rq, Ji C J 2 if and only if D 0 [Ji) D Dq[J 2 ]. 

Proof. It’s clear that J\ C J 2 implies D 0 [Ji] D Dq[J 2 \. Suppose J\ (f J 2 , so 
there is polynomial f £ J\ with f d J 2 . By Proposition 12.71 there is q £ Dq\J 2 ] 
with q(f) ^ 0, so Dq[J 2 \ <f_ D 0 [Ji\. ~ □ 

An immediate consequence of this corollary is that an ideal J C i?o is 
uniquely determined by its dual space Dq[J], 

Corollary 2.9. The dual space Dq[J] is homogeneous (respectively, monomial) 
iff the ideal J C Rq is homogeneous (respectively, monomial), i.e., generated 
by homogeneous elements with respect to filtration {m fc }fc>o (respectively, by 
monomials). 

Proof. Given a homogeneous (respectively, monomial) dual space L = Dq[J } of 
an ideal J C Rq it is straightforward to write down homogeneous (respectively, 
monomial) I C R such that Dq[I] = L. Namely, its homogeneous part of order 
k is the set of polynomials orthogonal to L k /L k ~ x \ for the monomial case, it is 
particularly explicit: a monomial x a belongs to I iff d a L. The extension IRq 
is determined by L uniquely according to Proposition 12.81 hence, IRq = J. □ 

Remark 2.10. One could easily extend the definition of homogeneous and 
monomial ideals to the local ring R y for an arbitrary point y £ C^: in par¬ 
ticular, an ideal is called monomial if it is generated by elements of the form 
(x - y) a , a £ (Z> 0 ) N . 

Macaulay dual bases allow for testing ideal membership at a solution [18] as 
stated in the following proposition. This can be readily generalized for homo¬ 
geneous ideals using the following corollary, Remark 12.41 and Proposition 12.81 

Corollary 2.11 (Lemma 11 of [3). A polynomial f £ R of degree d is a member 
of a homogeneous ideal I C R iff f is annihilated by Dq [/]. 

Proof. It follows from the proof of Corollary 12.91 that Dq[I] determines J = 
IR/m d+1 . Now, / £ / iff its image / £ J iff / is annihilated by Dq[I]. □ 

The statement of Corollary 12.111 corrects that of Theorem 4.6 of [15] where 
the assumption of homogeneity was missed as shown in [7]. The local mem¬ 
bership test without the assumption of homogeneity is a much harder task, 
addressed in m- 

2.2 Action of differentiation on the dual space 

An alternative characterization of the dual space can be given via Proposi¬ 
tion 12.131 There is a natural action of Rq on Dq by pre-multiplication. Specif¬ 
ically for q £ Dq and g £ Rq let g ■ q £ Dq denote the functional defined by 
(ff ' ?)(/) = 9 ( 9 /)- It can be checked that this gives Dq an R 0 -module struc¬ 
ture. The action of each variable Xi can also be considered as differentiating 
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functionals in Do by dt (up to normalization). Let a Xi : Do —> Dq denote the 
map defined by the action of Xi. 

Xi '■ Do —> Dq 

d a ^ d a ~ e \ (i = l,...,N), 

where d@ is taken to be 0 when any entry of /3 is less than zero. 

The following statements (from Proposition 12.121 to Corollary 12. 161) appear, 
perhaps in alternative phrasing, in many works addressing the duality at hand 
(see, for example, ns)- We collect the essential pieces, stated in our language, 
and complete with our own short proofs to guide reader’s intuition for this 
paper. 

Proposition 2.12. For a subspace L C Do the following are equivalent: 

• L is the dual space of some ideal Jl C Rq ■ 

• L is closed under differentiation by each variable: Xi ■ L C L for all 
1 < i < N. 

• L is an Rq- submodule of Do- 
Proof. For any L C Do define 

Jl = {/ S R 0 : q(f) = 0 for all q £ L}. 

If L is closed under differentiation, then Jl is closed under multiplication 
by each Xj, and therefore under multiplication by all monomials in R 0 . Express 
any g £ R 0 as g = c a x a . Then if / G J L and q £ L, q(gf) = c a q(x a f) 
and each term is zero, so gf £ Jl- Therefore Jl is an ideal and Dq[Jl\ = L. 
Conversely if L = D 0 [Jl] and q £ L then (x, • q)(f) = q{xtf) = 0 for all f £ Jl 
so Xi ■ q £ L. □ 

Consider the map 

Dual : {ideals of R 0 } —> {l?o“ su bmodules of D 0 } 

defined by Dual(J) = Do[J]- By Corollary 12.81 and Proposition 12.121 this map 
is a bijection. This provides another way to characterize the dual space. 

Proposition 2.13. For ideal J = (/i,..., /„) C Rq, let L be the maximal Ro~ 
submodule of D 0 that satisfies q{ff) = 0 for all q £ L and all 0 < i < n. Then 
L = D 0 [J]. 

Proof. Do[J] is closed under differentiation and satisfies q(fi) = 0 for all q £ 
Do[J] and 0 < i < n, so Dq[J] C L. The ideal Jl contains {/i,..., /„}, so 
J C J L which implies L C Dq[J]. □ 

Remark 2.14. For an ideal J £ Rq, the dual space Dq[J\ is finitely-generated 
as an Ro-module only when it is a finite dimensional vector space. If Dq[J] is 
generated by a single functional p, then J is exactly the apolar ideal of p (see, 
for instance, m for the definition). 
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A result of Proposition ^. 1 31 is that for I = (/i,..., /„}, a dual element q is in 
D 0 [I] if and only if q{fi) = 0 and Xj-q £ D 0 [J} for each 0 < i < n and 0 < j < N. 
Note that this leads to a completion scheme for computing D y [I] proposed in 
|19] , assuming y is in the vanishing set of I: 

D y[I] «- span c (<9°) 
for * = 1 —y k do 

D l[I\ {q e Dy I Xj ■ q £ D l ~ l [I] for all j = 1,... ,7V and q(fi) = 
0 for all i = 1,..., n} 

end for 

Moreover, the above algorithm makes apparent that if D y [/] = D * +1 [/] for some 
i > 0 then D l y [I\ is equal to all higher truncations, and so is equal to D y [I\. 
This gives an effective stopping criterion for computing D y [I\ when it is finite 
dimensional. 

Proposition 2.15. For ideals J\,Ji C Rq, 

• Do[Ji + J 2 ) = Dq[Ji\ H Do[J2\. 

• d 0 [Ji n J 2 ] = d 0 [Ji\ + Dq[j 2 \. 

Proof. The first statement follows from the definition of the dual space, as does 
Do[Ji] + Do[J 2 ] C Dq[Ji n J 2 ]. 

Let L = Dq[Ji] + Do[J 2 ]. It’s clear that L is an i?o-submodule, so it is 
the dual space of an ideal Jr.- The fact that Do[Ji] C L implies Jl C Ji and 
similarly Jl C Ji- Therefore -D 0 [Ji n J 2 ] C L. □ 

Corollary 2.16. If J\ and J 2 are homogeneous ideals of Rq, then the equality 
holds for the truncated dual spaces: 

-Dq[J r n J 2 ] = D$[J 1 ] + D d 0 [J 2 ], for all d £ N. 

Proof. This follows from the fact that if q £ D 0 [J] for a homogeneous J, then 
qdi the part of q of order d, is also in D 0 [J], □ 

Remark 2.17. For truncated dual space, in general, only one inclusion holds: 

Dg[Jin J 2 ] D D^Ji] + D*{J 2 \. 

However, because Dq[J\ n J 2 \ is finite dimensional, it follows that 

D*[J 1 nJ a \cD l 0 [J 1 ]+D l 0 [J 2 ] 

for l large enough. 

Example 2.18. Let I\ = {x\) and I 2 = (xi — x\) in R = C[xi,x 2 ]. Then 
D l[h] + D llh\ = span{l,d 2 }, 

Dl[h n/ 2 ] = span{l,9r,5 2 }, 

D o [ J i ] + D o [ j 2 ] = span{ 1, d 2 , d\ , d\ + d x } . 

There are strict inclusions 

Dhlh] + Dllh] C Dllh n I 2 ] c Dq [Jr] + D 2 0 [I 2 }. 
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2.3 Dual space of quotient ideals 

Recall that for g G R 0 , the map a g : D 0 D 0 denotes the action of g on D 0 
by pre-multiplication, or equivalently by “differentiation” with respect to g. 

Proposition 2.19. For all non-zero g G Rq, the map cr g : Dq — > Dq is surjective 
and kercr g = Dq [(<?)]. 

Proof. Note g ■ D 0 is closed under differentiation. If a g is not surjective, then 
g ■ Dq is the dual space of some non-trivial ideal / C Ho by Proposition 12.121 
Choose some non-zero / G I. Since gf^ 0, there exists some functional q with 
q(gf) ^ 0. Then g ■ q(f) = q{gf) ^ 0, which is a contradiction since g ■ q should 
annihilate /. To show ker a g = D 0 [(g)], if q G H 0 [(ff)] then g ■ q(f) = q(gf) = 0 
for all f G Rq. The only functional that is zero on all elements of Rq is the zero 
functional so g ■ q = 0. Conversely if q ^ D 0 [(g)} then g ■ q(f) = q(gf ) ^ 0 for 
some / G R, so g ■ q ^ 0. □ 

Theorem 2.20. Dq[I : (g)] = g ■ Dq[I], 

Proof. For / homogeneous, the statement is shown in [51 Theorem 22]. Here we 
consider the general case. 

If p G D 0 [I], then g ■ p(f) = p(gf) = 0 for all / such that gf G I. These are 
precisely the polynomials / in / : (g), and so g ■ p G Dq[I : ( 5 )]. 

For any q G D 0 [I : (< 7 )], because a g is surjective we can choose some p G D 0 
such that g ■ p = q. Then for all f G I : (g), we have q(f) = p(gf) = 0, so 

P e D 0 [g(I : < ff »] = D 0 [I n (g)} = D 0 [I} + D 0 [(g)}. 

Therefore p = p' + u for some p' G ZDo [T] and u G Do [(< 7 }]- Then q = g ■ p = 
g ■ p' + g ■ u but g ■ u = 0 so q G g ■ Dq [/] . □ 

3 Local Hilbert function and its regularity 

The Hilbert function of an ideal I G Rq provides combinatorial information 
about I that can be computed numerically using truncated dual spaces. 

3.1 Primal and dual monomial order 

Let > be a local monomial order (1 is the largest monomial), which we shall 
refer to as a primal order. For g = Yl a a aX a , a nonzero polynomial, the initial 
term with respect to > is the largest monomial with respect to > that has a 
nonzero coefficient, namely 

in>(< 7 ) = max]/ | a a ^ 0 }. 

For an ideal /, the initial terms of I with respect to > is the set of initial terms 
with respect to > of all the elements of /, namely 

in>(7) = {in>(/) | f G I}. 
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A monomial is called a standard monomial of I with respect to > if it is not a 
member of in>(7). 

We shall order the monomial differential functionals via the dual order: 

d a hdP & x a < x' 3 , 


the order opposite to >. 

The initial term in y(q) of q is the largest monomial differential functional 
that has a nonzero coefficient. The initial support of a dual space with respect 
to y is the set of initial terms with respect to >z of all the elements in the dual 
space (which can be considered as a subset of (Z> 0 ) N ). 

A dual basis that has distinct initial terms is called a reduced dual basis. 
Using a (possibly infinite dimensional) Gaussian elimination procedure, it is 
easy to see that any dual basis can be transformed into a reduced dual basis. 

Theorem 3.1 (Theorem 3.1 of jl6]). Let Iq be a O-dimensional ideal of Rq. 
The initial support of the dual space Dq [To] is the set of standard monomials for 
I = Iq D R, i.e., 


in b (A>[J 0 ]) = MA)[/]) = {d a | ar“ t in>(7)}. ( 2 ) 

Proof. Note that Dq[I] is finite dimensional. Choose a monic reduced basis B 
for Dq[I] such that the lead term of each element does not occur in any other 
element (using Gaussian elimination). 

Suppose d a € in^(Z? 0 [7]) so some p £ B has in>-(p) = d a . For any monic 
polynomial / with in>(/) = x a , f and p have no terms with the same exponent 
except their respective lead terms, so p(f) = 1 and / ^ I. 

Suppose d a in^(Do[/]). Let {pi,... ,p s } C B be the basis elements with 
d a in their monomial support. For each pi let Ci be the coefficient of d a and let 
< 9 / 3 * 

= m>(pi). The following polynomial 


f = x° 


£ 

i= 1 


has p(f) = 0 for all p £ B, and in >(/) = x a . By Proposition 12.71 / G Iq so 
x a G in>(/ 0 ) = in>(/). □ 

Corollary 3.2. For an ideal I C Ro, dime (Rq/(I + m fc+1 )) = dime -Dq [/] 
where m = ( X \,..., Xn) . 

Proof. Choosing a graded primal order >, a vector space basis for the quotient 
Ro/{I + m fe+1 ) is the set of monomials 

{x a | x a £ in>(J), and \a\ < fc}. 

By Theorem l3.1l this is corresponds to a basis for in^(I?Q [/]) which has the same 
dimension as Dj) [/]. □ 

We can extend Theorem 13.II to ideals of arbitrary dimensions. 
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Theorem 3.3. For an ideal I C R the monomial lattice is a disjoint union 
of in^ Dq[I] and in> 7. 

Proof. By Theorem 13.11 \ in>- Dq[I + m fc+1 ] = in>(7 + m A;+1 ), so then 

\ (J in h D 0 [I + m fc+1 ] = f| in>(7 + m fc+1 ). 

k k 

By definition lj fc in>- Dq[I + m fc+1 ] = in^ Dq [7], while by Lemma EU 
P|in>(/ + m fc+1 ) = in>(77 ? 0 H R) = in>(7). 

k 


□ 


3.2 Hilbert function and regularity index 

Definition 3.4. For an ideal I C Rq define the Hilbert function as 

f I + tn fe \ 

77/(7) = dim c (gr(7?o/t)fe) = dim c I ^ fc+1 1 

= dime (7?o/(7 + m fc+1 )) - dime (7? 0 /(7 + m fe )) . 

This is the same as 775'/j 0 // m , the Hilbert-Samuel function of the T?o-module 
Rq/I where Rq is filtered by {m fc }. 

The Hilbert function is determined by the initial ideal with respect to the 
primal monomial order (that respects the degree). 

Proposition 3.5. For an ideal I C Rq 

77/(7) = 77/ o (7) = t7 in >(/ n //)(7), for all k £ N. 

Alternatively, truncated dual spaces determine the Hilbert function. By 
Corollary 13.21 it can be seen that 

77/(7) = dimcT)o[7] — dime F 1 ^ 1 [7], for 7 > 0, 

where dime is taken to be 0 . 

For some m > 0 the Hilbert function is a polynomial in 7 for all 7 > m 
(see, e.g., B Lemma 5.5.1]), the Hilbert polynomial HP/(7). If the dimension 
of 7 C Rq is d , then HP/(7) is a polynomial of degree d — 1. In particular if 7 
is 0-dimensional then HP/(7) = 0 since Rq/I is finite dimensional. 

Definition 3.6. The regularity index of the Hilbert function is 

Po(I) = min{ m : Hj(k) = HP/(7) for all k > m}. 

The regularity index of an ideal is used as a stopping criterion for many 
algorithms which work iteratively by degree. In particular we will make use of 
it in Algorithm 15.11 
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4 Quotient ideals and eliminating dual spaces 

Section [2731 described the relationship between the dual space of an ideal Dq[I\, 
and the dual space of the quotient ideal by a principal ideal Dq[I : (c/)]. For 
applications (such as in Section[5]) it is useful to compute information even about 
the simplest case, where g = X\. We would like to find bases for the truncated 
dual spaces Dq[I : (xi)] but this proves difficult. 

Let > be a graded primal order on the monomials of the local ring Rq , and 
>~ be the dual order for the dual monomials of Dq. For any p £ Do, we must 
have ordin^(xi ■ p) < ordin^(p) — 1, since differentiation reduces the degree of 
each monomial by 1, but may also annihilate the lead term. Therefore taking 
the derivative of the dual space truncated at degree d+1 we have xi • Dq +1 [/] C 
Dq[I : (xi)]. Equality may not hold since there may be some functionals q £ 
Dq[I : (xi)] with q = xi ■ p for some p £ Dq[I] with lead term having degree 
higher than d + 1 and is annihilated by xi. In general, finding Dq[I : (xi)] 
from the truncated dual space of I may require calculating Dq[I] up to a high 
degree c. 

To overcome the difficulty of computing truncated dual spaces of colon ideals, 
we consider other filtrations of Dq corresponding to gradings on Rq other than 
the total degree grading. For A C {xi,..., x n } define ord .4 d a = Y^x-eA a *; the 
total order of all <9,; with Xi £ A. For general q £ Do define ord^ q to be the 
maximum order of the terms of q. 

Definition 4.1. Fixing A C {xi,..., xjv}, the eliminating truncated dual 
spaces of I are 

E$[I,A\={q£D 0 [I}:ord A q<d} 

for all d £ N. 

We often drop the word truncated when talking about eliminating dual 
spaces. 

The truncated dual spaces D d [I] give a filtration of D 0 [I] corresponding to 
the maximal ideal m of Rq 


Dq[I} = Dq[I + m d+1 ]. 

Similarly, the eliminating truncated dual spaces for A correspond to the ideal 
(A) in that 

E$[I,A]=Do[I+(A) d+1 ]. 

To see this, note that Eq[I,A] is the intersection of Dq[I\ with EofD,A\ = 
.Do[(A}“ +1 ]. By Pronositon 12.151 the intersection of these two dual spaces is 
D 0 [I + (A) d+1 }. 

For which ever grading of Rq (and corresponding filtration of Dq) is chosen, 
it is useful to pick a local order > (and corresponding dual order that is 
compatible with the grading. An order is compatible if for x“ £ (Ro)i and 
x^ £ (Ro)j with i < j then x a > x& . In the case of the total degree grading, 
such an order is a graded order. For the grading given by (A) a compatible 
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local order is an elimination order , eliminating the variables in A. In particular 
this is a block order in which the most significant block is a degree order on the 
variables in A and the second block is an arbitrary order on the variables not 
in A. Such an order ensures that p £ Eq [0, A] if and only if in>- p £ [0, A], 

Remark 4.2. Dual spaces offer analogs to many operations in elimination the¬ 
ory. The dual space of I n C[x m +i ,..., aqv] is equal to -Do[/]|a 1= o,..., 9 m =o- The 
eliminating dual E$[I,A\ is the dual space of I + (A), the variety of which is 
the intersection o/V(J) with the coordinate subspace in which the variables in 
A are zero. For a more detailed discussion in the case if homogeneous ideals see 

w- 

Let > be an local elimination order for aq,..., x m with dual order ^ and 
consider ring extension R' := C.(x m +i ,..., aqy)[aq, ■ ■ ■, x m ] A R- In this exten¬ 
sion with the corresponding dual order, the monomials in my Dq[IR'] are 
the monomials of in>- Dq[I] considering only the aq,... ,x m parts. These can be 
computed from in>- E$[I, {aq ,..., x m }} for sufficiently large d. 

Note that the eliminating dual space generalizes the usual truncated dual 
space since Eq[I,{xi, ... ,Xn}] = -Dq[/]. For general A, we have E$[I,A] D 
Dq[I\. Unlike Uq[/], the eliminating truncated dual space can be infinite¬ 
dimensional. 

Proposition 4.3. If I is an m-dimensional ideal that is in general position 
with respect to x i,.. •, x m then dime Eq[I, {x±, ..., x m }] < oo. 

Proof. For I satisfying these hypotheses the intersection of Y(I) with the space 
V(aq,..., Xm) is 0-dimensional. By Theorem 12.21 I + (x i,..., x m ) d+1 has dual 
space of finite dimension. □ 

In particular, if I is a curve, after a generic change of coordinates one can 
finitely compute its eliminating dual spaces for A = {aq}. The following propo¬ 
sition provides a method to compute eliminating dual spaces of the quotient 
ideal I : (xf) as well. 

Proposition 4.4. E$[I : (x\), {aq}] = x\ ■ E d+1 [I , {aq}] for all d £ N. 

Proof. Let ^ be a dual order on Dq eliminating x±. For any functional p £ Dq, 
either in> (p) is divisible by d\ , or p has no terms divisible by d\. In the first 
case, in^(a;i • p) = in^ (p) jd\. In the second case x\ ■ p = 0. Therefore, in the 
view of Theorem l2.201 any non-zero q £ E$[I : (aq), {mi}] must be the derivative 
of some p £ Eq +1 [I, {aq}]. □ 

This proposition is used in Algorithm 15.11 see Example 15.41 
Proposition 14.41 for curves does not hold in general (only a weaker Proposi¬ 
tion [475] does) and we are unable to use the eliminating dual spaces outside the 
specialized Algorithm 1-5.11 
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Proposition 4.5. 


Eq[I : (xi,... ,x m ),{xi,. .. ,x m }} D ^ Xj • £^o +1 [/, {xi,...,x m }] (3) 

2 = 1 

/or all d £ N. 

Proof. The inclusion ([3]) holds, since I : (x±,... ,x m ) = I '■ ( x i) and, by 
Theorem 12.201 

m m 

Dq[I : (xr,... , x m )] =^2d 0 [I : (Xj)] = ^2xf Dq[I], 

2=1 2=1 


□ 

Remark 4.6. Assuming it is finite, a basis for E$[I, {xi,..., x m }] can be com¬ 
puted by finding a basis of the dual space of I + (xi,..., x m ) d+1 . The dual 
space of a 0-dimensional ideal can be efficiently computed for example with the 
algorithm of m or others. 


5 Detecting embedded points on curves 

The general problem of detecting an embedded component can be formulated 
as follows: 

Consider an ideal I C R and a prime ideal Po 3 I. Let P\,...,P r D 
I be associated ideals of I such that \fP[ C P 0 . 

Given generators of I together with general points y o € V(P) and 
y, £ V(Pj) [i = 1, ■ ■ • , r) determine whether P is an associated prime 
of R/I. 

The problem that we solve here is more special: we consider the case when 
the variety is locally a curve, namely, dim yo I = 1. That means dimP^ = 1 
for i / 0 and V(Po) = {z/o} is a point that may or may not be an embedded 
component. 

Let an ideal I be given by its generators F and suppose, without a loss of 
generality, that the point in question is yo = 0. Let the 1-dimensional primary 
components in the problem be Pi,..., P r with Vi = V(Pj) containing the origin. 
Saturating I by the ideal (xi) eliminates all the components of I that contain 
(xi). After a generic linear change of coordinates, we may assume that no V r is 
contained in the hyperplane X\ = 0 except for Vq = V(Qo)> so I : (xi)°° ^ I if 
and only if the origin is an embedded component. 

This leads to the following algorithm that employs the eliminating dual 
spaces. 

Algorithm 5.1. B = IsOriginEmbeddedlnCurve(J) 
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Require: I, a 1-dimensional ideal of R in regular position relative to X\. 
Ensure: B = “origin is an embedded component of I”, a Boolean value. 

r <- Po{i); 

™ <- MoCO; 

k <r- max(r, m — 1); 

E<^E*[I,{x i}]; 

return X\ ■ E C Eg '[Id* i}] 

Here /x 0 (I) denotes the multiplicity (or degree) of / at the origin. For I 
a curve, note the (local) Hilbert polynomial of / is the constant polynomial 
HP/(fc) = no(I). To compute po(I) and po(I) we refer to the algorithm given 
in m which produces the Hilbert function of I from a set of generators, and 
in the process the Hilbert regularity index and the Hilbert polynomial of I. 

The following two lemmas are used in the proof of correctness of Algo¬ 
rithm 15.11 

Lemma 5.2. Suppose ideals I,Jc Rq satisfy I C J and dim qJ/I is finite. 
Then I = J if and only if 


-Do” 1 ! 7 ] = -DS _1 [ J ] 


where r = max{/?o( 7 ), Po(J)}- 

Proof. Since I C J, to show I = J it is enough to show that Hj(k) = Hj(k) 
for all k > 0. Because dime J/I is finite, HP/ = HPj so the Hilbert functions 
agree for k > r. If additionally Dg~ l [I] = Z?g -1 [J], then the Hilbert functions 
also agree for 0 < k < r. □ 

Lemma 5.3. If J is a one-dimensional monomial ideal that is saturated at the 
origin (J = J : m°°), then 


Po{J) < Po(J) ~ 1- 

Proof. We consider a monomial cone decomposition of the standard monomials 
of J. For monomial m £ Rq and a set of variables v = {x ^,..., Xi k } the 
monomial cone C m:V is 

C m , v ■= { x...x a Pm\ (ax,..., a k ) £ N fe }. 

A monomial cone decomposition of Rq/J is a finite list of pairs 

(mx,vx),..., (m s ,v s ) 

such that the standard monomials of J are a disjoint union of the cones 

fArai ,vi i * • ■ j f-Oris ,v a • 

The dimension of a cone C mj1) is defined to be the size of v. A cone decompo¬ 
sition is closely related to the Hilbert function of J: the maximum dimension 
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of a cone in the decomposition is the dimension of the ideal; the number of 
maximal dimensional cones is the multiplicity po(J); and the maximum de¬ 
gree of the monomials mi,... ,m s bounds the regularity po(J). For J a one¬ 
dimensional monomial ideal saturated at the origin, there is a cone decomposi¬ 
tion (mi, {x^}),..., (m s , {xi a }) of Rq/J consisting only of dimension 1 cones. 

Modify this decomposition slightly by letting to' := rrij\ Xi .—\, the monomial 
obtained from rrij by removing Xi j . The cones C m ' Ul ,..., C m > iVa also have the 
standard monomials of J as their union, but are generally not disjoint. To 
prove the proposition, it is sufficient show that for all d > po(J) — 1 each cone 
contains exactly one monomial of degree d and these monomials are distinct, 
and therefore Hj(d) = po (J). 

Let Mfe := {m' | ij = k}. Note that J2k l-^fcl = Po(J)- For each k , Mk is 
closed under differentiation. This follows from the fact that Mk is the set of 
standard monomials of 7r( J) where tt : C[xi,..., xjv] C[xi,..., Xk, ■ ■ ., xjv] is 
the projection sending Xk to 1. If Mk has a monomial of degree d , it also has 
at least one monomial of each degree < d. Therefore 

max degm < \Mk\ — 1 < po(J) — 1. 

mGMfc 

So for d > fiQ (J) — 1, each cone contains a monomial of degree d. 

Suppose two cones in the decomposition intersect, so n = m' x“ = mjx^ for 
some j ^ l and Xi j ^ x,,. Then x“ divides mj so a < degmj. 

deg n = deg m' + a < deg m' + deg m[ < \M it | + \M ij \ - 2 < p 0 - 2. 

No two cones have a monomial in common of degree d > po(J) — 1. □ 

Proof of correctness of Algorithm 1 5. il By Proposition 14.41 xi • E = F?q - 1 [/ : 
(xi),{xi}]. If this dual space is not equal to Eq~ 1 [I,{x i}] then / : (xi) ^ I. 
This implies there is an embedded component at the origin. 

Suppose instead Xi • E = Eq^ 1 [I,{xi}]. We will use Lemma Hn2l to prove 
that / : (xi) = I. The truncated dual space of degree r — 1 is contained in the 
eliminating dual space of degree r, so Dg _1 [/ : (xi)] = Hg _1 [/]. We know that 
I C I : (xi). Because they differ by at most a zero-dimensional component, 
dimc(/ : (xi))// is finite. 

Finally it must be shown that k > max(p 0 (d), po(I : (xi))). It is clear that 
k > p 0 (I)- To show k > po(I : (xi)), let J = in(/) : m°°, which has the same 
Hilbert polynomial as I and I : (x\) and satisfies 

ffi(-I) Q in (I : (xi)) C J, 

Hi > H I: ( Xl ) > Hj. 

By Lemma [5731 po(J) < Po(I) ~ 1- Since U/./ xl \ is sandwiched between Hj and 
Hj , once they stabilize to po(I), so must This implies the regularity of 

I : (xi) is bounded by k. □ 
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Example 5.4. Let I = { x 2 — z 3 ,y — z 2 ) C C[x,y,z] which defines a curve 
in C 3 with a singular point at the origin. The deflation algorithm from m 
will identify the origin as a possible embedded component. Note that po(I) = 
1 , po (I) = 2 and no irreducible component o/V(/) is contained in the plane 
x = 0. To test whether the origin is embedded, we compute the eliminating dual 
El[I , {x}]. This is the set of all dual functionals with all terms having d x -degree 
< 1 . 

Eq[I, {x}] = spanjl, d 2 + d y , d z , d x ,d x d 2 z + d x d y ,d x d z }, 
x-Eq[I,{x}] = span{l ,d 2 + d y ,d z }. 

Since x- Eq[I, {x}] = £g[J, {x}] we conclude that the origin is not an embedded 
component of I. 

Example 5.5. For this example we compute with an implementation of Al¬ 
gorithm \5J] in Macaulay2. Let I be the ideal of the cyclic4 system, generated 
by 

{Xi + X 2 + X 3 + X 4 , X’iX 2 + X 2 X 3 + X3X4 + X4X1, 

X2X3X4 + X1X3X4 + XiX 2 X4 + XiX 2 X 3 , XiX 2 X 3 X4 — 1 }. 

I is a curve with several singular points, which are discovered using the algorithm 
described in m up to some numerical precision. One such point is p = 

(- 0.0000000000000000122 + 1 . 0000000000000002221 , 

-0.0000000000000001128 + 0.999999999999999889 i, 

0.0000000000000000459 - 0.9999999999999998891, 

-0.0000000000000000935 - 1.0000000000000000001), 

approximately (1,1, —1, —1). Let I' denote the ideal obtained from I by a random 
affine change of coordinates that fixes p. This ensures that I' is in general 
position with respect of x\. Using the algorithm of ms, the regularity index is 
pp(I') = 2 and the multiplicity is p P (I') = 1, so k = max(2, 0) = 2. 

Computing E 3 [/',{xi}] and x\ ■ E 2 [I',{x i}] the dimensions are 3 and 2 
respectively, so they are not equal. Therefore the point being approximated by p 
is an embedded component of I. The code for this example can be found at ffjf l. 
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